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ABSTRACT 

Current ramp designs used in RORO operations have been determined to be 
structurally inadequate in Sea-State 3. The main reason for this 1s that when the ramp is 
loaded with two vehicles, the relative motions between the ship and the discharge facility 
induce stresses above yield. The objective of this thesis is to formulate the problem and 
present results that will enable the design of proper isolators. The problem is formulated 
in terms of the hydrodynamic interaction between adjacent bodies along with structural 
coupling. The applicability and limitations of strip theory approaches are established 
through comparisons with three dimensional hydrodynamic analysis data. An analytic 
model of the ship-ramp-barge system 1s developed and tested. The results indicate the 
validity of the approach and establish a procedure that may be utilized for the design of 


passively controlled isolators. 
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I, INTRODUCTION 

Over the years, there has been very little research conducted on vessels with 
objects in the vicinity. Of the research conducted, most of it was centered on the analysis 
of supply ships mooring alongside oi] platforms. With the increased emphasis on littoral 
warfare, it is becoming more important to be able to discharge supplies from ships to 
shore without a readily available port facility. Therefore, it is desired to develop supply 
ships with the ability to offload cargo to a roll-on, roll-off discharge facility (RRDF) in 
the open ocean. 

To accomplish this, a ramp 1s mounted from the ship to the RRDF. However, due 
to hydrodynamic forces on the ship and barges, the mounting connections are exceeding 
their design limitations. Therefore, it has become necessary to develop a method to 
minimize these forces on the pinned connections. To accomplish this, the hydrodynamic 
forces on the ship and barge must be accurately modeled. Then, these hydrodynamic 
forces must be adapted to determine the magnitude of the forces that act on the pinned 


connection so that a passive damping system may be produced in the future. 
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Il. BACKGROUND 


A. MOTIVATION 

In recent years, the United States military has taken an increasing interest in 
developing the ability to offload equipment and machinery from roll-on, roll-off (RORO) 
vessels in the open seas. Operational requirements for equipment transfer are set at sea 
state three while maintaining a total dynamic acceleration of no more that 25% of the 
static gravitational constant. In order to accomplish equipment transfer, a connected ramp 
system has been established between the RORO and a transport barge. The operational 
requirements stated above should hold in the case of two tanks located in the middle 
point of the ramp. This corresponds to an envisioned worse case scenario of one tank 
breaking down in the middle of the ramp and the second tank providing assistance. This 
coupled problem of multiple rigid bodies combines elements of hydrodynamic 
interactions and structural response. The philosophy of a design and analysis strategy is 
how to ensure proper connectivity between the different bodies such that cargo motions 
and ramp stress levels are minimized. 


B. OBJECTIVE 


The objective of this thesis is to formulate and present preliminary results for the 
problem of reducing ramp stress levels via reductions in ramp motion inputs. This will be 
done by coupling the hydrodynamic analysis of the adjacent floating bodies to a 
structural model of the ramp. Passive/semi-active isolator design will then be used as the 


basis for the synthesis of the entire model in the frequency domain. 


C. BACKGROUND 


A great deal of the analysis of hydrodynamic study of multiple body systems 
revolves around the wave effects upon vertical cylinders in waves. Most of this research 
surrounds the desire to accurately predict the affects of waves on different parts of a 
structure to ensure geometric stability. Much of this research is motivated to produce 
continued stability in offshore 01] producing stations during drilling operations. 

However, there has not been significant analysis performed when two independent bodies 


are in close proximity to one another. [Ref. 9: p. 725] 


Developing the equations of motion of a solitary vessel in various incident wave 
conditions is fairly routine. It simply involves defining the added mass characteristics 
and associated stiffness and damping matrices, which the ship will possess due to its 
design. However, the effects that waves will have on ships in close proximity to one 
another are very different. This is because the diffracted and radiated waves and their 


respective hydrodynamic forces will affect the adjacent vessels. 


The equations of motion of vessels in close proximity to one another are derived 
in much the same manner as independent vessels. However, when multiple vessels are 
close to one another, additional terms must be incorporated into the equations of motion 
to ensure validity. The importance of these additional terms depends on relative size of 
each vessel, how each vessel is positioned with respect to the other, and overall distances 


between them. [Ref. 9] Taking these criteria into account, an additional added mass term 


and an additional damping term must be utilized to accurately predict the hydrodynamic 


response of each vessel. 


Research has indicated that there are multiple techniques that could, in principle, 
predict the numerical values of the additional terms necessary to accurately predict the 
hydrodynamic motion of a vessel. One of these techniques is the use of what is known as 
the macroelement technique. This method approximates each of the vessel’s cross- 
sectional areas with a step curve. The flow field around each area is then subdivided into 
macroelements and approximated. This technique uses the method of separation of 
variables and Laplace’s differential equations for the flow field. Once the velocity 
potentials for the diffraction and radiation features for each element are established, a 
Fourier series is developed that meets the free surface and ship body boundary 
conditions. To obtain the remaining coefficients, the Galerkin method is utilized. This 
method is effective as 1t does not produce many numerical irregularities and saves in 
computer processing time. [Ref. 1] Due to its approximations, this method may be 
gees for certain cases and it suffers from the basic limitations of a two-dimensional 


analysis, most notably end effects and the existence of standing waves. 


Another technique that is widely studied is that of finite element approximations. 
In this approach, the area around the vessel is modeled using boundary dampers/infinite 
elements or eigenfunction expansion/boundary integral techniques. [Ref 5: p 358] Using 
standard assumptions that wave and body motions are small, the complex velocity 
potential is obtained for the diffracted and radiated wave potentials. Using data regarding 


complex velocity potentials; information about linearized free surface, bottom, radiation, 


and body boundary conditions; and hydrodynamic forces and coefficients, a final 
equation is developed. This equation integrates and sums the added mass, damping, and 
stiffness matrices of each vessel and its effects on the other vessel. Furthermore, it solves 
for the forces on each vessel due to waves and any associated hydrodynamic effects from 
the other vessel. This method is effective for various wave conditions as long as the 
problem assumptions are met. However, proper nodal placement and meshing is 
essential for accurate results. Due to the requirement for very fine meshing in some 


problems, computer resources may restrict the determination of accurate results. [Ref 5] 


Techniques based on strip theory formulations are easier to solve than the 
previous methods while producing an accurate wave force and linked body prediction. In 
this approach, there are two floating vessels attached alongside one another. The 
equations of motion are generated using known data for each vessel and coupled together. 
When coupled, each equation of motion from one vessel effects that of the other. Vessel 
interaction is approximated by assuming that as a wave strikes a vessel, the wave is 
reflected an infinite number of times between vessels. This produces two additional 
equations for the wave interaction forces on each vessel. To determine the wave 
interaction forces, strip theory, which assumes that changes in the longitudinal gradients 
of the flow field around a vessel are much smaller than corresponding changes in the 
transverse direction, is utilized. The wave forces are a combination of the Froude-Krylov 
force and the diffraction wave force. Additional equations for added mass, damping, and 
stiffness for each vessel are needed to represent the interaction effects. Taking the 


equations generated by wave interaction between the vessels and those equations for 


added mass, damping, and stiffness and placing them back into the original equation of 
motions for each vessel, a new equation of motion is developed which will depict the 

motion of the vessel to various wave loading. This technique has shown to be effective 
in all sea directions as long as the assumptions of linearized motions and slender bodies 


are accurate. [Ref 2] 


Another related technique that is commonly used to solve for the additional terms 
is that of the boundary integral method. This method first divides the vessels into 
sections. The radiation velocity potential function is determined using Green’s theorem 
and the Fredholm integral equation of the second kind, which is then applied to each of 
the divided sections. Once the radiation velocity potential function is developed, an 
interpolating shape function can be formed. Finally, a set of linear equations are 
developed which produce a set of algebraic equations for each of the radiation velocity 
potential equations throughout the vessels. This method produced accurate results in all 


types of wave loading directions at varying frequencies. [Ref 3] 


D. CONTENTS 


The boundary integral technique is the method that most researches utilize when 
obtaining the numerical solution to the added mass, damping, and stiffness terms in the 
equation of motion of a vessel. For this reason, it was adopted as the method of choice for 
this study. Chapter III presents the theoretical foundation for the hydrodynamic equations 
of motion. This is followed by the development of the algebraic constraint equations that 
were incorporated in the original differential equations of motion in order to model the 


structural coupling between the bodies. Modeling of the connecting ramp 1s also 


7 


included in Chapter III, followed by the results of this work in Chapter IV. Finally, the 
conclusions from this research and some recommendations for follow-up work are 


summarized in Chapters V and VI. 


Ul. PROBLEM FORMULATION 
A. INTRODUCTION 
When examining the connection between the ship and the barge, the points where 
the ramp attached to the barge and the ship exhibit induced forces on the pinned 
connection. The equations to develop these forces are formulated by considering heave, 
pitch, and roll forces upon the ship and the barge. 


B. HYDRODYNAMIC EQUATIONS OF MOTION 


To obtain the force equations, it was necessary to define a standard coordinate 
system. The system used was the same one utilized in the reference 14, Principles of 
Naval Architecture, in which n1=surge; noa=sway; n3=heave; n4=roll; ns=pitch; and 
Ns=yaw. The three principle motions that affect the problem were heave, roll, and pitch. 
Relative surge, sway, and yaw motions were considered negligible and were not analyzed 


further. 


The other assumption that was made at this point was that the only appreciable 
affect upon ship movement would be in the vertical direction. Therefore, the equation for 


absolute vertical motion of a vessel at a given point was found to be 


E,=n3+,y- xn, (1) 
Where y and x is the motion at any point on the body in the x and y planes. &; is the 


complex local amplitude at x, and y in heave. 


After defining the governing equation that would control the motion of each 


vessel, a free-body diagram was produced of the ramp. This diagram assumed that there 


9 


are three spring connections mounting the ramp to the ship and barge. Each connection 
produces a displacement at the pinned joint defined as &), &), and 3. The first € 
represents the connection for the ship and the last two &’s are for the barge connections. 
The next assumption was that each spring would exert a force at the connection as a 
force-couple system. This force-couple system represents a single force at each pinned 


location. Therefore, the equation given for the force at each location is defined as 


F =G.é. 
us (2) 


where the subscripts 1 and j reference the particular force-couple system at a given 
location. G is the matrix of influence coefficients, which is a transfer function. From (2), 


the following matrix was obtained: 
F, ee Titel 
Ga Gs Galnes (3) 
F, Gen Gesell. 

Using this process, it was then necessary to find the equations of motion for 
heave, roll, and pitch. This was done using (4) where A, B, and C are the added mass, 
damping, and stiffness matrices for the ship and barge. The j terms are for the 6 
directions of motion in which only heave, roll, and pitch were considered. F; and F;” are 
the complex amplitude of the exciting force components due to incident and diffracted 


waves, and F; are the connection forces at the ramp, to be determined. 


6 
sy e aa, A,)+ io,B,+C, he =Fl+F °+F (4) 
k=) 


Solving these equations for the three relevant motions to this problem, the final equations 


of motion were found to be: 


10 


Asm, + Ass, = F, in Heave (5) 
Aun, = F, in Roll Cop 


Assn, + Ass. = F, in Pitch (G7) 


The various terms 1n these equations are defined as: 


n,(t)=7,e'"' => Heave (8) 
n,(t)=7,e'*" => Roll (9) 
n(t)=7,e'° => Pitch (10) 


Each of the internal terms in the equations of motion are all matrices. These matrices are 
characterized as A, B and C and represent the added mass of the vessel, the damping 
coefficient, and the stiffness. This implies that each vessel’s hydrodynamic properties are 
modeled as a spring-mass-damper system. The first number in the following equation’s 
subscripts refers to the direction (heave, pitch, or roll) and the second number is with 
respect to a particular direction. For example, A35 in (12) would be the heave added mass 


due to pitch. 


Ass =-@, (A+A,,) +i,B,, +C,, (ian) 
Ass = 1 [Ax + As} +10.B,, + Cy (12) 
Ass =—00, (I, +A.) +i@,B,, +C,, “3 ) 
Ass =~ [We +4) +1@2B, +G, (14) 
M=-olier A, | tiaB, +C,, o>) 


The terms in the added mass equations are the specific equations of mass, damping, and 
stiffness for the particular vessel. Due to inherent symmetrical relations within the 


vessel, the following additional assumptions are utilized: 


I 


Cy, = C's (16 ) 
BGs > @oe. (L7®) 


Ay = A 5; (18) 


The force terms in the equations of motion are simplified to make individual force 
calculations easier. This leads to the following equations that produce the resultant 


heave, roll, and pitch forces: 


F,=F, + F, + F,’ (19) 
Bae ae es (20 ) 
Fa Sa, eye a ied (21) 


The F3', F4', and F;' terms are the heave, roll, and pitch forces that act on the ship or barge 


from the motion of the ramp. 


The equations of motion in (5) through (7) are body independent. Consequently, 
they are used to describe the motions of the ship and barge. Therefore, a subscript ‘S’ or 
‘B’ is used to distinguish the motions of each vessel. Furthermore, equations (5) through 
(7) fail to recognize that as two objects are in close proximity to one another, there will 
exist a radiative effect between the two vessels. This results from waves striking one 
vessel, bouncing off the hull and then striking the vessel alongside. These reflected 
waves effect heave, pitch and roll of the adjacent vessel. Combining the vessel 


distinguishing subscripts and the radiative matrices yielded the following equations. 


Heave, ship > A338 fag + Ass,s ie + Bss.a ae + Bss.3 oe + Base Was = Fy; (223) 


Pitch, ship => As3.s 7755 + 455.5755 + Bs3.8155 + Bss.ets_ + Bsieny, = Fes 3) 


te 


(24 ) 


Roll, ship => Ass as a Bases 5 + Bus.eNs 5 5 5 Bas. 5 = Fy 


Heave, barge > A338 1.5 + As5,8 17s 5 +Bustys + Bass 1s.s + Bast gs — F 3 ( 2) ) 
Pitch, barge => As. 7 + Ass.B Tews + Bsus T)5,s + Bss.s 5.5 + Bss.s tas =F» ( 26 ) 
Roll, barge => ows Tis Beare fe Bs ie Bias iets = Fy, C2.) 


The B terms are the radiative matrices between the barge and ship while the subscript S 


and B designate the ship or barge. The radiative matrices are defined with the following 


equations: 
Baas =-0'°D,,, +iwE,,, <2) 
Bus =-@°D,y,, t+i@E, 5 (29) 
Bae = OUD». 5 4 Ona. (30 ) 
ae = -a° D4 5 To eae (31) 
Bass = -0° Dy, tioky , (32) 
ee =-0°D,,,+i@k 4 5 (35) 
Bs.e = -o° Ds , +iwk,, , (34) 
Bssp = -0° D5, t+iWE sy 5 (5) 
Bss.p = —@7 Dos p +iWE ss p (36) 


D and E are the added mass and damping influence coefficients due to the radiated waves 
from one vessel to another. Similar equations for the radiative matrices are obtained for 
the ship. Due to symmetry Ex =E, and Dyf=D,j. An example of one of these 
symmetric equations 1s the following: 


D552 = Ds3 5 (37) 


IS 


Using a free-body diagram, (5) through (7) were expanded allowing the 
hydrostatic and ramp forces on each vessel to be separated. The terms Fy‘, F2‘, and F;' 


were defined as f}, fo, and f3 respectfully. The following revised equations of motion are: 


Heave, ship > A3351]3.5 + Ass ee + B3357]3.5 + Bss.8 745 + Bass ee = ee qh (38 ) 
Pitch,ship=> As3s1)5.5 + Asss Ts 5 + Bssp ies + Bsas Nas + Bsa), 5 = ae —fix, (39) 
Roll, ship=> Aass74,s + Bas 7.5 + Bass ore + Bas es = i + hy (40) 
Heave, barge > A3.81)3,5 + Ass. ie + Bast}; + Bass ns + B3ssT)s 5 = ee ae face (41 ) 
Roll, barge=> Ass. Fs + Bas). + Bassas + Bass ]s. = Bes eel ( 4° ) 


Pitch, barge => A538 1] 5 + Ass.81]s 5 + Bs3,s ae + Bss.s7]s5 + BsssT)s.s = Be — f,X,—f3X; (43 ) 


This allowed for the effects on the ship and barge at the pinned connections due to the 
direct forces from the ramp as well as moment forces due to ramp movement. The 


superscript H is to show the force acting on the hull of the respective vessel. 


Once the above hydrodynamic equations of motion are solved, the values for f}, 
f>, and f; can be established. However, all of the n terms in the above equations are 
unknown and all the f terms are also unknown. Therefore, the values for the n terms are 
determined by defining additional equations. 7 is thus defined as a function of f), f2, and 


f; for the ship and barge as: 
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lo = fo (44) 


Nas = Ha fi + Ha (45) 
age elles ec (46 ) 
(EO fo OL eS Oe (47) 
a5 =D fee Die Ee 5, (48) 
i.» =m, +ooe fo, i) 


From these equations, it is possible to solve the six linear equations for each respective u 
and v in terms of various known added mass matrices, known physical positions on the 
vessel, and the forces acting on the pinned connections as well as the hydrodynamic wave 


forces. 


Once each u and v term is solved, every n is considered a known value in terms 
of f;, f2, and f3. Each n term has a contribution to the total force acting where the ramp 
attaches to the barge or the ship. Therefore, the total force acting on a pinned joint 
includes forces from heave, roll, and pitch. This produces three combined displacement 


equations for each connection that are: 


C1 See es) eee i iG 0) 


oo i ee ne = Koy fy thax fy + Kr (51) 
$3 = Tae 1 sa V3 = rales = ky fy + Ka; fs + K 30 (52) 


It was stated previously that €), 2, and €3 are the combined displacements that act 


at each of the pinned connections. Now the equations are defined by relating the forces 


1D 


acting on the pinned connections from the above k terms using the known formulas for n 
and equating portions of the € equations that correspond to the respective x or y terms. 


Thus, the k terms become: 


ky = Hay + May ~ Hs) Rig = [oe ee Lea (53 ) 
Ky9 = 39 + Vg Yq — V59X3 Koz = 33 + Uy Vo — U53X2 Kg = V3q + Vy V2 — Vs9X, (34) 
Kyy = Vy +0 V3 — 52%; Kyy = Us; + Ug, V3 — Us3X; R59 = Vyq + VyV3 — VsoX3 (5353 


These k terms are then set equal to the transfer functions that relate the forces at one 
pinned connection to the other forces from each pinned connection, and how the various 


forces interact through the ramp system. This leads to the equations: 


Kiet Hid) = Oo Ol. 1en. (56 ) 
Kop + ka fo t+ hu hs = Gan fit Gn fr + Guhs (57) 
Ky t ky fr t ky fs = Gun fi t+ Gun fy, + Ga; fs (58 ) 


Rearranging equations (56) through (58) and placing into matrix for yields: 


k, =O -Gi, ~G I —kyo 
-Gy, ky -Gy ky -G3 | fF. |=] — ho ( 59 ) 
—GF k;, -G;, k,; -G;; ts —kyo 


This first matrix is labeled as matrix [G,]. Solving for the forces yields: 


ti — ki 
es Ic, J — k a9 ( 60 ) 
I 3 — Kk 5 


The three forces are easily solved using Matlab or any other engineering computational 


program given the numerical values of the added mass, damping, and stiffness matrices. 
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C. RAMP EQUATION OF MOTION 


1. Introduction 

To accurately model the ship barge combination, it is important to incorporate the 
ramp connecting the two vessels because there will be interaction effects between the two 
vessels transmitted through the ramp. The ramp is modeled here as a rigid body, 
although the procedure is general enough so that it can be generalized to include 


flexibility effects. 


2. Derivation of the Equation of Motion 

To obtain the equations of motion for the ramp, it was first necessary to examine 
the ramp as a Separate entity. First, the ramp is a three-degree of freedom system in the 
heave, pitch, and roll directions. The ramp is assumed to be connected to the ship and 
barge at three locations. See appendix A for figure of LMSR/barge connection. Each of 
these locations has a spring connection and the overall ramp has a mass of M. It was 
further assumed that the ramp is a homogenous device with uniform properties 
throughout. Using these assumptions, the equations of motion were derived using the 
method of summing kinetic and potential energies of the ramp and then using Lagrange’s 


Equations to derive the equations of motion. 
First, it is necessary to define Lagrange’s equations. The Lagrangian 1s defined as 


b= (61 ) 


where L is the Lagrangian, T is kinetic energy, and U is the potential energy of the ramp. 


The Lagrangian is then used to find the value of the forces at each degree of freedom. 
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The equation for this is written below as equation (62). This equation takes the partial 
derivative of the Lagrangian with respect to the velocity of one of the degrees of freedom. 
The partial derivative of the Lagrangian with respect to position is then taken and 
subtracted from the time derivative of the first value. This produces the force at one of 


the degrees of freedom. 


£{ 2). 2h. (62) 


After the Lagrangian is defined, the kinetic and potential energies for the ramp are 
found. Looking at the ramp from a view approximately 45 degrees above the horizontal, 
it was assumed that there are three pinned connections. These are at the two corners that 
would attach the ramp to the barge and the third connection is from a kingpin in the 
center of the ramp that attaches the ramp to the ship. These connections produce a force 
upward due to the associated springs. Visualizing the ramp in some deformed state, it is 
possible to define the change in position due to heave, pitch, and roll. The heave is 
simply the distance y upward from the original axis. The pitch is defined as a positive 
upward angle from the new horizontal axis on the ship side of the ramp. Roll is the 


rotation about the x-axis from the new pitch angled ramp. Kinetic energy is defined as 


e 2 2 2 . 2 4 
T= .5M y +.5M6 4.5My (63 ) 


In this equation, M is the mass of the ramp and y, 8, and PV are the three individual 


displacements which act on the ramp throughout the ramp motion. 
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Potential energy is defined as the stored energy of the springs throughout the 
system. Damping of the ramp was not included in the determination of the potential 
energy as proportional damping was considered to be in effect. Therefore, the value of 
the damping coefficients will be proportional to the spring coefficients by some set 
amount. This allows easier calculation of the potential energy. In the potential energy 
equation, U is the potential energy, k is the spring coefficient, A is the length of the ramp, 


and B is the width of the ramp. Y, 0, and are the same as previously defined. 


2 


2 2 
A, A”. Bee a. B 
U =.5k,| y+—sin © | +.5k —-—sinO+—sin ¥ | +.5k ~——sin © —- —sin ¥ 
(v 2 ; 0 2 2 0 iid i (64) 


CA 


Once the kinetic and potential energies are defined, the determination of the 
equations of motion of the system is fairly straightforward. Set up the Lagrangian, take 


derivatives, and equate terms and place in matrix form. The force equations were found 


to be 
M = bf y+ Asin | - ka{y- Sin © + Ssin v|- ba[ y- Asin © - sin | =a (65 ) 
em. @k ie in® ee ech One Brey By i ee ier = 
Go aD 2 Pe Cee weed OG?) 
M ¥ —~-cos wk,{ y- Asin Open wv )- Zo whs{ » Oe amt Aa, Ou 
2 2 2 2 2 2 
It should be noted that there is not a term for ’ in the k; part of the equation. This is 


because as the ramp rolls, the changes in the ¥ angle for the spring on the centerline is 


zero. This also causes the Fw force to not have a k; term as the derivative does not have a 


YY term in it. 
A A A B B A A A B B 
; © —<C, -C, on OF —G +-C, +—-C, —C, +—C, _ —K, -K, -K; — K, +— K, +— K, — K, +— K, 
EF MOO y ° - 2 5 2 z y ee 2 2 ; 2 2. y 
2 a A A A A A A AB. AB ee A A A A A A AB AB 
nh M 0 e —C, +—C, +—C, —C, — CC, — C, —C, ae O — K, +—K, +— K, — K,— kK, —— K; —— K, —K, © 
Fyl {0 0 M\y 2 2 | 2 |F 
BB AB, AB B B B AB AB BB 
= a nog ee 3 ha tf, —K; re yk ee 


(68 ) 
Equation (68) above shows the end result matrix for the three directional forces on the 
ramp. This matrix 1s produced after deriving the equations of motion for the ramp, 
adding the proportional damping values and linearizing the all three equations using the 
small angle approximation. This allowed sin@ and sin'¥ to be approximated as simply @ 


and ‘Y and cos@ and cos'¥ to be approximated as 1. 


Once the final matrix for the ramp equations of motion was derived, these forces 
representing the position of the ramp in terms of y, 9, and ‘¥ had to be correlated to be in 
terms of F;, F2, and F3, the forces at each of the connecting locations and their associated 
changes in absolute vertical motion. First, the forces on the ramp were written in terms 


of F;, F2, and F3. This produced the following matrix where all symbols are as previously 


Fas | 1 | F, 

P = > = = fe 

Fy 0 oo = "Bet (69 ) 
2 


stated. 


The next step was to state the changes in position of the ramp connecting locations in 


terms of y, 9, and 


Hite SD 
a ae |e 
*3 j fe a 
Z 2 


In this matrix, X;, X2, and x3 are the vertical displacements for each of the pinned 
connections. The final relation must be to have Fy, Fe, and Fy in terms of only the 


displacements. The resulting equation is as follows: 


F, y 
F, |=[-o?[M ]+ jof[c]+[K]] © 
Fy y 
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(71) 


In this matrix, the M, C, and K matrices are the mass, damping, and spring stiffness 
matrices defined earlier. Frequency is defined as m and is the frequency of wave impact 


from the hydrodynamic information. 


Once the relations are made between the pinned connections and the equations of 
motion, a relation can be established to relate the displacements of the connections to the 


ramp equations of m © 


A 
1 = 0 
x, Ss a Y oy l 1a 
- - -A 
=e ——— IE Poa i aaa 
|i eaten eee 2 Rm eeeeee ers 
2 2 Z 2 


Known values are the length and width of the ramp, the M, C, and K matrices, and the 
hydrodynamic forces on the pinned connection. Placing all this information together 
results in the final matrix describing the change in vertical displacement of the pinned 


connection. 
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IV. RESULTS 


A. LARGE MEDIUM-SPEED ROLL-ON ROLL-OFF VESSEL ANALYSIS 


1. Introduction 

The comparison between the strip theory and full three-dimensional (WAMIT) 
program solvers are presented in Figures 1 through 48. Each of the values in the figures 
has been non-dimensionalized in a consistent manner displayed in Table 1. Furthermore, 
all pertinent ship parameters are shown in Table 3. All symbols used are identified in 


Table 2. 











Ship/Barge Coefficient WAMIT 
Non-dimensionalization 
A 


o 


B53, Pitch to Heave Damping Bs3/opL” 
C33, Heave Stiffness C33;*pgA,/pgL 
F3, Heave Exciting Force F;/pgAL 
F5, Pitch Exciting Moment F./pgAL 
Phase Angle 


Table 1. Table of Various Non-Dimensionalizations 







2. Data Analysis 


Figures 1 through 28 display the added mass and damping results for the large 
medium speed roll-on roll-off vessel (LMSR) data. The first 8 figures portray data 
obtained from the strip theory solver and WAMIT when over 2000 panels were used in 
the 3-D analysis. The next 8 figures show the same added mass and damping curves 


comparing strip theory and WAMIT but in these curves, there were only about 1000 
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panels for the ship and barges. Figures 21 through 28 show the effects due to wave 
radiation effects from the barge on the LMSR. Overall, the results show that relatively 
good agreement is obtained between the two computer solving techniques, especially in 
the heave mode. However, oscillatory behavior is observed for the WAMIT program at 
high frequencies (low periods), especially when fewer panels are utilized. This is due to 
discretization effects in the solver. To overcome these problems, a relatively high 
number of panels are required to accurately represent the hydrodynamic interactions 
between the LMSR and the small Roll-on, Roll-off Discharge Facility (RRDF). The strip 
theory solver does not take such interactions into consideration. Therefore, it is not 


susceptible to high frequency “chattering”, unlike the 3-D solver. 


DEFINITION 


Mass of the ship or barge in pounds. 


SYMBOL 








M 


L Length of the ship in feet 


Force due to gravity in pounds-force per second-squared. 


Waterplane area produced by the strip theory solver in square feet. 
Wave amplitude. 
Frequency of wave impact upon the ship or barge in Hertz. 


Ship’s displacement in cubic feet. 


SS | Density of Seawater in pounds-mass per cubic feet. Considered to be 1.9904 Ib/ft’. 


Table 2. List of Applicable Symbols and Their Definitions 








As mentioned, the 3-D solver can overcome the high frequency oscillations by 
utilizing a more refined panel scheme. However, the very large number of panels that is 
required by the 3-D solver incurs a significant cost in both computational time (more than 


16 hours for the 2000 panel case on a high-end CPU) as well as the time needed to 
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prepare the data file. By contrast, the strip theory solver results were obtained after only 
a few seconds of execution time on an older generation Pentium CPU. Of course, the 
limitations of the strip theory solver are that the interaction effects are not capable of 
being modeled. By comparing the strip theory and WAMIT solvers, it was clear that 
modeling of the added mass and damping characteristics of the main ship could be 
achieved both accurately and efficiently through a combination of the two-dimensional 
and three-dimensional hydrodynamic solvers. 

The results for the hydrodynamic wave excitation forces and phase angles for the 
LMSR are presented in Figures 25 through 60. These wave forces range from zero 
degrees (head seas) through 315 degrees at 45-degree intervals. This provides complete 
360-degree coverage for wave force analysis for the LMSR. In these plots, overall 
agreement between the strip theory solver and WAMIT is outstanding. As will be noted 
in more detail when the barge data is analyzed, there are some slight differences between 
the port and starboard wave directions in the 3-D results. These deviations are due to 


interactions between the barge and the LMSR. 






[Dimension _—~«| SSR 


Ramp Location (Feet From FP) X=301.2, Y=71.3 39.5 (From Hinged Side Edge) 
8.8(From Centerline of Causeway 
to Port) 


Table 3. Applicable Ship, Barge, and Ramp Dimensions 

















2D 


3. Figures 


aor Heave Added Mass - LMSR 





5.5 [ 


A33 


3.5 | 


\ Pa ¥ | 





0 2 4 6 8 10 HZ 14 16 18 
Period (sec) 


Figure 1. Plot of Strip Theory and Wamit 
Program Comparisons for Heave Added Mass 


for the LMSR. 
a0. Heave Damping - LMSR 

5 ee 
4.5} > | 
4° A 
35° a 
, ’ | 

3t . 

B VA 
< 2.5 ¥ | 
2° | — > strip | 4 

——_wamit 

1.5f x 4 








0 Z 4 6 8 10 12 14 16 18 
Period (sec) 


Figure 2. Plot of Strip Theory and Wamit Program 
Comparisons for Heave Damping for the LMSR. 


26 


ean Pitch Added Mass - LMSR 
3.5 TOO OS ll Lr 


- 


: | 
3 ; | 


| 


0 2 4 6 8 10 12 14 16 18 
Period (sec) 





Figure 3. Plot of Strip Theory and Wamit Program 
Comparisons for Pitch Added Mass for the LMSR. 
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Figure 6. Plot of Strip Theory and Wamit Program 
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Figure 8. Plot of Strip theory and Wamit Program 
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Figure 10. Plot of Strip Theory and Wamit Program 
Comparisons for Heave Damping for the Old LMSR 
Data. 
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Figure 11. Plot of Strip Theory and Wamit Program 
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Figure 12. Plot of Strip Theory and Wamit Program 
Comparisons for Pitch Damping for the Old LMSR 
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Figure 14. Plot of Strip Theory and Wamit Program 
Comparisons for Heave/Pitch Damping for the Old 
LMSR Data. 
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Figure 15. Plot of Strip Theory and Wamit Program 
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Figure 16. Plot of Strip Theory and Wamit Program 
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Figure 17. Plot of Multiple Barge Comparison for 
Heave to Pitch Added Mass. 
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Figure 18. Plot of Multiple Barge Comparison for 
Heave to Pitch Damping. 
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Figure 19. Plot of Multiple Barge Comparison for 
Pitch to Heave Added Mass. 
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Figure 20. Plot of Multiple Barge Comparison for 
Pitch to Heave Damping. 


35 








18 


18 


x1 o Heave Added Mass of LMSR due to Heave Added Mass of Barge 








2 | We | 
4 
2) ; 1 
1.5) = —— 1st Barge onLMSR q 1 
Pe —=- 2nd Barge on LMSR a | 
1 - 1 
05° ' ee . 
| an = 8 
O} =" | Esp ae 
-0.5+ as 
=f = 
=f Oar 
! a | 
2 3 | 
2 “ 6 8 10 12 14 16 18 
Period(in seconds) 
Figure 21. Plot of Multiple Barge Comparison for 
Heave Added Mass of LMSR due to Heave Added 
Mass of Barge. 
x 10° Heave Damping of LMSR due to Heave Damping of Barge 
4 
| ‘i. o£ 
3+ J 
3 a 
QL — 
! oD, —=- 1st Barge onLMSR 
| ge? —2-— 2nd Barge onLMSR 
a a 
© 1} ¢ fa 5 


Oo 
3 
Se 


a y 
aq 
| | 
<2 a ee a he ee re | 
2 4 6 8 10 12 14 16 18 
Penod(in seconds) 


Figure 22. Plot of Multiple Barge Comparison for 
Heave Damping of LMSR due to Heave Damping of 
Barge. 
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Figure 23. Plot of Multiple Barge Comparison for 
Pitch Added Mass of LMSR due to Pitch Added Mass 
of Barge. 
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Figure 24. Plot of Multiple Barge Comparison for 
Pitch Damping of LMSR due to Pitch Damping of 
Barge. 
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Figure 25. Plot of Multiple Barge Comparison for 
Heave Added Mass of LMSR due to Pitch Added 
Mass of Barge. 
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Figure 26. Plot of Multiple Barge Comparison for 
Heave Damping of LMSR due to Pitch Damping of 
Barge. 
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x 40 Pitch Added Mass of LMSR due to Heave Added Mass of Barge 
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Figure 27. Plot of Multiple Barge Comparison for 
Pitch Added Mass of LMSR due to Heave Added 
Mass of Barge. 
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Figure 28. Plot of Multiple Barge Comparison for 


Pitch Damping of LMSR due to Heave Damping of 
Barge. 
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Figure 29. Plot of Strip Theory and Wamit Program 
Comparisons for Heave Exciting Force at an Incident 
Wave Angle of 0 Degrees for the LMSR. 
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Figure 30. Plot of Strip Theory and Wamit Program 
Comparisons for Heave Exciting Force Phase Angle 
at an Incident Wave Angle of 0 Degrees for the 
LMSR. 
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Figure 31. Plot of Strip Theory and Wamit Program 
Comparisons for Pitch Exciting Force at an Incident 


Wave Angle of 0 Degrees for the LMSR 
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Figure 32. Plot of Strip Theory and Wamit Program 


Comparisons for Pitch Exciting Force Phase Angle at 
an Incident Wave Angle of 0 Degrees for the LMSR. 
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Figure 33. Plot of Strip Theory and Wamit Program 
Comparisons for Heave Exciting Force at an Incident 
Wave Angle of 45 Degrees for the LMSR 
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Figure 34. Plot of Strip Theory and Wamit Program 
Comparisons for Heave Exciting Force Phase Angle at 
an Incident Wave Angle of 45 Degrees for the LMSR. 
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Figure 35. Plot of Stnp Theory and Wamit Program 
Comparisons for Pitch Exciting Force at an Incident 
Wave Angle of 45 Degrees for the LMSR 
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Figure 36. Plot of Strip Theory and Wamit Program 
Comparisons for Pitch Exciting Force Phase Angle at 
an Incident Wave Angle of 45 Degrees for the LMSR. 
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Figure 37. Plot of Strip Theory and Wamit Program 
Comparisons for Heave Exciting Force at an Incident 
Wave Angle of 90 Degrees for the LMSR. 
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Figure 38. Plot of Strip Theory and Wamit Program 
Comparisons for Heave Exciting Force Phase Angle 
at an Incident Wave Angle of 90 Degrees for the 
LMSR. 
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Figure 39. Plot of Strip Theory and Wamit Program 


Comparisons for Pitch Exciting Force at an Incident 
Wave Angle of 90 Degrees for the LMSR. 
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Figure 40. Plot of Strip Theory and Wamit Program 
Comparisons for Pitch Exciting Force Phase Angle at 
an Incident Wave Angle of 90 Degrees for the LMSR. 
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Figure 41. Plot of Strip Theory and Wamit Program 
Comparisons for Heave Exciting Force at an Incident 
Wave Angle of 135 Degrees for the LMSR. 
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Figure 42. Plot of Strip Theory and Wamit Program 


Comparisons for Heave Exciting Force Phase Angle at 
an Incident Wave Angle of 135 Degrees for the LMSR. 
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Figure 43. Plot of Stip Theory and Wamit Program 
Comparisons for Pitch Exciting Force at an Incident 
Wave Angle of 135 Degrees for the LMSR. 
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Figure 44. Plot of Strip Theory and Wamit Program 


Comparisons for Pitch Exciting Force Phase Angle at 
an Incident Wave Angle of 135 Degrees for the LMSR. 
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Figure 45. Plot of Strip Theory and Wamit Program 
Comparisons for Heave Exciting Force at an Incident 
Wave Angle of 180 Degrees for the LMSR. 
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Figure 46. Plot of Strip Theory and Wamit Program 


Comparisons for Heave Exciting Force Phase Angle at 
an Incident Wave Angle of 180 Degrees for the LMSR. 
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Plot of Wave Exciting Force for Pitch at 180 Degrees 
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Figure 47. Plot of Strip Theory and Wamit Program 


Comparisons for Pitch Exciting Force at an Incident 
Wave Angle of 180 Degrees for the LMSR. 
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Figure 48. Plot of Strip Theory and Wamit Program 
Comparisons for Pitch Exciting Force Phase Angle at 
an Incident Wave Angle of 180 Degrees for the 
LMSR. 
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Figure 49. Plot of Strip Theory and Wamit Program 
Comparisons for Heave Exciting Force at an Incident 
Wave Angle of 225 Degrees for the LMSR. 
Plot of Wave Exciting Force Phase Angle for Heave at 225 Degrees 
200 
150 
100 + 
50 - 
ea) 
50 | 
-100 
150} 





Penod(in sec) 


Figure 50. Plot of Strip Theory and Wamit Program 
Comparisons for Heave Exciting Force Phase Angle at 
an Incident Wave Angle of 225 Degrees for the 
LMSR. 
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Figure 51. Plot of Strip Theory and Wamit Program 
Comparisons for Pitch Exciting Force at an Incident 
Wave Angle of 225 Degrees for the LMSR. 
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Figure 52. Plot of Strip Theory and Wamit Program 
Comparisons for Pitch Exciting Force Phase Angle at 
an Incident Wave Angle of 225 Degrees for the 
LMSR. 


=a 


Amplitude 


Phase Angle(in Degrees) 


Plot of Wave Exciting Force for Heave at 270 Degreess 








0.07 ; 
: eee 
0.06 } ae J 
0.05 - * u 
A 
| 4° —— Stip- | 


Oo 
Oo 
W 
* 
eS ee eee 





ih 
0.02 j 
# 
f 
0.01 - a 4 
! fet 
1 
i —— 
o 2 4 6 8 0 2 4 6 418 


Penod(in sec) 


Figure 53. Plot of Strip Theory and Wamit Program 
Comparisons for Heave Exciting Force at an Incident 
Wave Angle of 270 Degrees for the LMSR. 
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Figure 54. Plot of Strip Theory and Wamit Program 
Comparisons for Heave Exciting Force Phase Angle at 
an Incident Wave Angle of 270 Degrees for the 
LMSR. 
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Figure 55. Plot of Strip Theory and Wamit Program 
Comparisons for Pitch Exciting Force Phase Angle at 
an Incident Wave Angle of 270 Degrees for the 
LMSR. 
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Figure 56. Plot of Strip Theory and Wamit Program 
Comparisons for Pitch Exciting Force Phase Angle at 
an Incident Wave Angle of 270 Degrees for the 
LMSR. 
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Figure 57. Plot of Stip Theory and Wamit Program 
Comparisons for Heave Exciting Force at an Incident 
Wave Angle of 315 Degrees for the LMSR 
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Figure 58. Plot of Strip Theory and Wamit Program 
Comparisons for Heave Exciting Force Phase Angle 
at an Incident Wave Angle of 315 Degrees for the 
LMSR. 


54 


Plot of Wave Exciting Force for Pitch at 315 Degrees 
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Figure 59. Plot of Strip Theory and Wamit Program 
Comparisons for Pitch Exciting Force at an Incident 
Wave Angle of 315 Degrees for the LMSR 
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Figure 60. Plot of Strip Theory and Wamit Program Comparisons for Pitch Exciting 
Force Phase Angle at an Incident Wave Angle of 315 Degrees for the LMSR. 
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B. ROLL-ON, ROLL-OFF DISCHARGE FACILITY (BARGE) ANALYSIS 

1. Introduction 

The hydrodynamic characteristics of the Roll-On, Roll-Off Discharge Facility 
(RRDF) are presented in Figures 61 through 106. These figures present a very detailed 
analysis of the barge and describe the effects of the large medium speed roll-on roll-off 


vessel (LMSR) on the barge. 


2. Data Analysis 


Figure 61 illustrates the added mass in heave for a barge unit by itself and also the 
leading (first) and trailing (second) barges located on the starboard side of the ship. The 
results displayed in all the figures are in nondimensional form as described in Table 1. 

It should be noted that the characteristic length which is utilized to nondimensionalize the 
values for the barges are that of the ship and not of the barge. The reasoning for this is to 
ensure all numbers are nondimensionalized using the same method, which allows for a 
more accurate comparison of results. 

The results displayed are those for the three-dimensional (3-D) solver only, which 
modeled the RRDF as a square platform with shallow draft. The actual dimensions for 
the RRDF are given in Table 3. Due to the lack of a precise draft, the results generated 
by the use of the strip theory solver were determined to be erroneous, exhibiting large 
divergences from the 3-D solver for many different values. Therefore, data was 
compared against a stationary barge labeled “1 Barge” in the figures. This generated the 
basis for establishing the relative significance of hydrodynamic proximity effects. The 


results of Figures 61 through 72 quantify the proximity effects on the heave and pitch 
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added mass and damping. These figures demonstrate the difference from the baseline 
data of the single barge to the barges adjacent to the LMSR is up to 30 percent. The 
irregular oscillatory behavior of the plots at lower periods is attributed to artificially 
generated high frequency interactions due to panel discretization effects. 

Results from the wave exciting forces are shown in Figures 73 through 106. 
These figures display the non-dimensional heaving force and pitching moment 
amplitudes and phase angles, for the headings ranging from Zero to 360 degrees at 
intervals of 45 degrees. Headings of zero degrees indicate waves coming directly from 
the bow of the RRDF, 90 degrees from the starboard beam, 180 degrees from the stern, 
and 270 degrees from the port beam. The results clearly indicated the significance, often 
two or more times higher, of the proximity affect on overall exciting forces. The very 


large effect of proximity forces was most clearly noted in the following figures: 


a. Proximity Induced Exciting Force 

Significant effect on pitch for beam seas as evidenced in Figure 83. Here the 
pitch force is essentially zero for the single barge case which is expected due to hull 
geometry. Therefore, any non-zero forces acting on the barges are entirely due to the 
presence of the LMSR. Similar results are found when the wave angle is from 270 


degrees in Figure 99. 


b. Sheltering/Magnification Effect 


It should be noted that both the heave and pitch exciting forces are significantly 
larger when the waves are from 90 degrees compared to when the waves are from 270 


degrees. The reason for this disparity is due to the incident waves being reflected off the 


a7 


starboard side of the LMSR and then impacting upon the RRDF’s which again create a 
magnification effect in the 90 degree case. By contrast, when the wave direction is from 
270 degrees the incident waves are scattered by the port side of the ship that produces a 


sheltering effect on the RRDF’s from wave action. 
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3. Figures 
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Figure 61. Plot of Single Barge and Multiple Barge 
Comparisons for Heave Added Mass. 
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Figure 62. Plot of One Barge and Multiple Barge 
Comparisons for Heave Damping. 
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Figure 63. Plot of One Barge and Multiple Barge 
Comparisons for Pitch Added Mass. 
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Figure 64. Plot of One Barge and Multiple Barge 
Comparisons for Pitch Damping. 
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x 10 Heave Added Mass of Barge due to Heave Added Mass of LMSR 


2 ‘Ss a 


15) Ee 
—— ist Barge p 
—3- 2nd Barge : 








2 4 6 8 10 12 14 16 18 
Period(in seconds) 
Figure 65. Plot of Multiple Barge Comparisons for 
Heave Added Mass of Barge Due to Heave Added 
Mass of LMSR. 
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Figure 66. Plot of Multiple Barge Comparisons for 


Heave Damping of Barge Due to Heave Damping of 
LMSR. 
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x 10” Pitch Added Mass of Barge due to Pitch Added Mass of LMSR 
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Figure 67. Plot of Multiple Barge Comparisons for 
Pitch Added Mass of Barge Due to Pitch Added Mass 
of LMSR. 
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Figure 68. Plot of Multiple Barge Comparisons for 
Pitch Damping of Barge Due to Pitch Damping of 


LMSR. 
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Figure 69. Plot of Multiple Barge Comparisons for 
Heave Added Mass of Barge Due to Pitch Added 
Mass of LMSR. 
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Figure 70. Plot of Multiple Barge Comparisons for 
Heave Damping of Barge Due to Pitch Damping of 
LMSR. 


63 


B113 
Oo 
- fe 
A 


x 10 Pitch Added Mass of Barge due to Heave Added Mass of LMSR 











() 
| | 
1 + 
Rh stag Sao 
| fT | fs | +— 1stBarge | 
i | 3 |—S- 2nd Barge | 
0.5} ' | : -| 
ho! i _ 
LAG at. a eee 
i | a oy - 
Oo; Ay A | | i 
| | L 
“1 | ee 
id, 
-0.5 Pd 7 | 
o 
Lr , | Ge 
: 
-1.5 
2 4 6 8 10 12 14 16 18 


Penod(in seconds) 


Figure 71. Plot of Multiple Barge Comparisons for 
Pitch Added Mass of Barge Due to Heave Added 
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Figure 72. Plot of Multiple Barge Comparisons for 
Pitch Damping of Barge Due to Heave Damping of 
LMSR. 
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Figure 73. Plot of One Barge and Multiple Barge 
Comparisons for Heave Exciting Force at an Incident 
Wave Angle of 0 Degrees. 
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Figure 74. Plot of One Barge and Multiple Barge 
Comparisons for Heave Exciting Force Phase Angle 
at an Incident Wave Angle of 0 Degrees. 
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Figure 75. Plot of One Barge and Multiple Barge 
Comparisons for Pitch Exciting Force at an Incident 
Wave Angle of 0 Degrees. 
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Figure 76. Plot of One Barge and Multiple Barge 
Comparisons for Pitch Exciting Force Phase Angle at 
an Incident Wave Angle of 0 Degrees. 
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Figure 77. Plot of One Barge and Multiple Barge 
Comparisons for Heave Exciting Force at an Incident 
Wave Angle of 45 Degrees. 
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Figure 78. Plot of One Barge and Multiple Barge 


Comparisons for Heave Exciting Force Phase Angle at 
an Incident Wave Angle of 45 Degrees. 
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Figure 79. Plot of One Barge and Multiple Barge 
Comparisons for Pitch Exciting Force at an Incident 





Wave Angle of 45 Degrees. 
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Figure 80. Plot of One Barge and Multiple Barge 
Comparisons for Pitch Exciting Force Phase Angle at 
an Incident Wave Angle of 45 Degrees. 
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Figure 81. Plot of One Barge and Multiple Barge 
Comparisons for Heave Exciting Force at an Incident 
Wave Angle of 90 Degrees. 
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Figure 82. Plot of One Barge and Multiple Barge 
Comparisons for Heave Exciting Force Phase Angle at 
an Incident Wave Angle of 90 Degrees. 
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Figure 83. Plot of One Barge and Multiple Barge 
Comparisons for Pitch Exciting Force at an Incident 
Wave Angle of 90 Degrees. 
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Figure 84. Plot of One Barge and Multiple Barge 
Comparisons for Pitch Exciting Force Phase Angle at 
an Incident Wave Angle of 90 Degrees. 
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Figure 85. Plot of Multiple Barge Comparisons for 
Heave Exciting Force at an Incident Wave Angle of 
135 Degrees. 
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Figure 86. Plot of Multiple Barge Comparisons for 
Heave Exciting Force Phase Angle at an Incident 
Wave Angle of 135 Degrees. 
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Figure 87. Plot of Multiple Barge Comparisons for 
Pitch Exciting Force at an Incident Wave Angle of 135 
Degrees. 


Plot of Wave Exciting Force Phase Angle for Pitch at 135 Degrees 


200 i - er | 








| 

150; | 

100 | | | 

—~— Barge 1, With LMSR 

sal —— Barge 2, With LMSR | 

| 

0 i + YN, } 

-50 | me eee. 

-100, | 

-150 | 

| | | 

-200 ; 

0 4 6 fio 812) 14 16 43 


Period(in sec) 


Figure 88. Plot of Multiple Barge Comparisons for 
Pitch Exciting Force Phase Angle at an Incident Wave 
Angle of 135 Degrees. 
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Figure 89. Plot of Multiple Barge Comparisons for 
Heave Exciting Force at an Incident Wave Angle of 
180 Degrees. 
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Figure 90. Plot of Multiple Barge Comparisons for 
Heave Exciting Force Phase Angle at an Incident 
Wave Angle of 180 Degrees. 
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Figure 91. Plot of Multiple Barge Comparisons for 


Pitch Exciting Force at an Incident Wave Angle of 180 
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Figure 92. Plot of Multiple Barge Comparisons for 
Pitch Exciting Force Phase Angle at an Incident Wave 


Angle of 180 Degrees. 
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Figure 93. Plot of Multiple Barge Comparisons for 
Heave Exciting Force at an Incident Wave Angle of 
22> Deorees, 


Plot of Wave Exciting Force Phase Angle for Heave at 225 Degrees 
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Figure 94. Plot of Multiple Barge Comparisons for 
Heave Exciting Force Phase Angle at an Incident 
Wave Angle of 225 Degrees. 
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Figure 95. Plot of Multiple Barge Comparisons for 
Pitch Exciting Force at an Incident Wave Angle of 
225 Degrees. 
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Figure 96. Plot of Multiple Barge Comparisons for 
Pitch Exciting Force Phase Angle at an Incident 
Wave Angle of 225 Degrees. 
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Figure 97. Plot of One Barge and Multiple Barge 
Comparisons for Heave Exciting Force at an Incident 
Wave Angle of 270 Degrees. 
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Figure 98. Plot of One Barge and Multiple Barge 
Comparisons for Heave Exciting Force Phase Angle 
at an Incident Wave Angle of 270 Degrees. 
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Figure 99. Plot of One Barge and Multiple Barge 
Comparisons for Pitch Exciting Force at an Incident 
Wave Angle of 270 Degrees. 
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Figure 100. Plot of One Barge and Multiple Barge 
Comparisons for Pitch Exciting Force Phase Angle at 
an Incident Wave Angle of 270 Degrees. 


78 


Phase Angle(in Degrees) 


Plot of Wave Exciting Force for Heave at 315 Degrees 





—— 1 Barge 


| —— Barge 1,WithLMSR | oo i 
9 —~ Barge 2, With LMSR_ ee 
Oa | 
4 ( 
os ) 
© : or 
s A cas - 
= 3 (a) oy 
E a - 
t A 
2t : 





16 


10 14 18 


Period(in sec) 


12 


Figure 101. Plot of One Barge and Multiple Barge 
Comparisons for Heave Exciting Force at an Incident 
Wave Angle of 315 Degrees. 
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Figure 102. Plot of One Barge and Multiple Barge 


Comparisons for Heave Exciting Force Phase Angle at 
an Incident Wave Angle of 315 Degrees. 
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Figure 103. Plot of One Barge and Multiple Barge 
Comparisons for Pitch Exciting Force at an Incident 


Wave Angle of 315 Degrees. 
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Figure 104. Plot of One Barge and Multiple Barge 
Comparisons for Pitch Exciting Force Phase Angle at 


an Incident Wave Angle of 315 Degrees. 
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Figure 105. Plot of One Barge and Multiple Barge 
Comparisons for Heave Exciting Force at an Incident 
Wave Angle of 90 Degrees and 270 Degrees. 
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Figure 106. Plot of One Barge and Multiple Barge 
Comparisons for Pitch Exciting Force at an Incident 
Wave Angle of 90 Degrees and 270 Degrees. 
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C. RAMP ANALYSIS 


1. Introduction 


The forces on the ramp exerted from the hydrodynamic motions of the LMSR and 
barge are presented in Figures 107 through 13la. These figures show the overall forces 
on the pinned locations that connect the ramp to the LMSR and the barge. To derive 
these forces, the data utilized was that obtained from WAMIT as this data produced 
excellent results when compared with actual wave tank model tests. The ramp is 
assumed to be a rigid body for analysis and stiffness of the spring connections is assumed 
to be proportional to the hydrostatic heave stiffness of the LMSR or barge, depending on 


which pinned connection is being analyzed. Table 4 shows the dimensions of the ramp. 


Ramp Length 
(in Feet, Pivot to Pivot) 


| ss Ramp Width (in Feet, Overall) == 


Ramp Weight 96.94 ee Tons 


Table 4. Ramp Dimensions 





Description 










2. Data Analysis 

Figures 107 through 122 show a plot of the ramp force magnitudes for the three- 
pinned connections. It is clear that the pinned connection that receives the largest peak 
magnitude force is the third king pin which connects the ramp to the barge. The 
connection that receives the largest force throughout the entire range is the first pinned 
connection. The third connection receives up to 5 times the magnitude of force as 
compared to the pinned connections attaching the ramp to the LMSR at the peak 


frequency. The first king pin has a force about 1.5 times greater than the third connection 
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throughout the entire range. The forces on the barge’s king pin are largest when the 
waves are from 90 and 270 degrees. 

When looking at the individual king pin for the LMSR and the barge, 
comparisons were made between the actual force at the pinned connection and the 
hydrodynamic heave force acting on the vessel. This data is displayed in Figures 123, 
125, and 127 where the actual force which each king pin experiences is higher than the 
hydrodynamic heave force on the vessel at higher wave frequencies. The force that the 
ramp connections experience drops off at longer periods. 

Figures 107 through 122 and 129 through 131a show the changes in overall 
excitation force on the pinned connection if different spring stiffnesses are assumed to 
exist. As can be seen in Figures 129, 130, and 131, the larger the value of the spring 
constant, the greater the magnitude of force which each king pin observes. This implies 
that a lower stiffness value for the spring connection is one method of minimizing 


excitation forces on the pinned connections. 
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Figure 107. Plot of Ramp Exciting Force on the First 
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Figure 108. Plot of Ramp Exciting Force on the 1“ 
Pinned Connection at Various Wave Angles with a 
Spring Constant of 0.1. 


84 


Magnitude of Force 


0.002 + 
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Figure 109. Plot of Ramp Exciting Force on the First 
Pinned Connections at Various Wave Angles with a 
Spring Constant of 1. 
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Figure 110. Plot of Ramp Exciting Force on the 
Second Pinned Connections at Various Wave Angles 
with a Spring Constant of 0.01. 
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Plot of Ramp Force Magnitudes at Various Angles for 2nd Pinned Connection, K=0.01C33 



























Magnitude of Force 





0.01 - ) ) 
0.009 r | | 
| : 
0.007 7 
| PEAR 
| a 90 Degrees 
0.006 /-—— 135 Degrees 
ee 180 Degrees | 
—o— 225 Degrees 
0.005 ® ~«— 270 Degrees | 
| ie | 
0.004 + ant, : 
1 / i ie 
0.003 + Mie % 3 
| ¢ / > — | 
0.002 | . a » - oe 1 
0.001 oS ae 
| : a a 
0 | 
2 8 10 12 14 16 18 





Penod (in Seconds) 


Figure 110a. Plot of Ramp Exciting Force on the 
Second Pinned Connections at Various Wave Angles 
with a Spring Constant of 0.01. 
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Figure 111. Plot of Ramp Exciting Force on the 2"° 
Pinned Connection at Various Wave Angles with a 
Spring Constant of 0.1. 
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Plot of Ramp Force Magnitudes at Various Angles for 2nd Pinned Connection, K=.1C33 
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Figure 111a. Plot of Ramp Exciting Force on the 2" 
Pinned Connection at Various Wave Angles with a 
Spring Constant of 0.1. 
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Figure 112. Plot of Ramp Exciting Force on the 
Second Pinned Connection at Various Wave Angles 
with a Spring Constant of 1. 
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Plot of Ramp Force Magnitudes at Vanous Angles for 2nd Pinned Connection, K=C33 
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Figure 112a. Plot of Ramp Exciting Force on the 
Second Pinned Connection at Various Wave Angles 
with a Spring Constant of 1. 


Plot of Ramp Force Magnitudes at Various Angles for 3rd Pinned Connection, K=0.01C33 
0.08 





| 
0.07 | ] 


0.06 - 
| | | 
® ! | 
: 0.05 + 1} 
u. | 
‘ 0.0 : 
® 0.04; ooo 
3 | —— 90 Degrees | | 
S —--— 135 Degrees 
«6 0.037; _-&- 180 Degrees a 
= | —=— 225 Degrees | 
| -— 270 Degrees | 
0.02 j J 








2 7 ae 8 10° «12 14 jou mis 
Period (in Seconds) 


Figure 113. Plot of Ramp Exciting Force on the 
Third Pinned Connection at Various Wave Angles 
with a Spring Constant of 0.01. 
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Plot of Ramp Force Magnitudes at Various Angles for 3rd Pinned Connection, K=0.01C33 
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Figure 113a. Plot of Ramp Exciting Force on the 
Third Pinned Connection at Various Wave Angles 
with a Spring Constant of 0.01. 
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Figure 114. Plot of Ramp Exciting Force on the 3" 
Pinned Connection at Various Wave Angles with a 
Spring Constant of 0.1. 


89 





Magnitude of Force 


Magnitude of Force 


Plot of Ramp Force Magnitudes at Vanous Angles for 3rd Pinned Connection, K=.1C33 


—+— 90 Degrees 
—— 135 Degrees 
—S- 180 Degrees 
| —+- 225 Degrees 
2/70 Degrees 








Penod (in Seconds) 


Figure 114a. Plot of Ramp Exciting Force on the 3 
Pinned Connection at Various Wave Angles with a 


Spring Constant of 0.1. 
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Figure 115. Plot of Ramp Exciting Force on the Third 


Pinned Connections at Various Wave Angles with a 


Spring Constant of 1. 
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Figure 115a. Plot of Ramp Exciting Force on the 
Third Pinned Connections at Various Wave Angles 
with a Spring Constant of 1. 
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Figure 116. Plot of Ramp Exciting Force on all Three 
Pinned Connections at a Wave Angle of 90 Degrees 
with a Spring Constant of 0.1. 
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Figure 116a. Plot of Ramp Exciting Force on all 
Three Pinned Connections at a Wave Angle of 90 
Degrees with a Spring Constant of 0.1. 
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Figure 117. Plot of Ramp Exciting Force Phase 
Angle on all Three Pinned Connections at a Wave 
Angle of 90 Degrees with a Spring Constant of 0.1. 
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Figure 118. Plot of Ramp Exciting Force on all 


Pinned Connections at a Wave Angle of 90 Degrees 
with all Spring Constants 0.1 of the LMSR Stiffness. 
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Figure 118a. Plot of Ramp Exciting Force on all 
Pinned Connections at a Wave Angle of 90 Degrees 
with all Spring Constants 0.1 of the LMSR Stiffness. 
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Figure 119. Plot of Ramp Exciting Force on all Three 
Pinned Connections at a Wave Angle of 135 Degrees 
with a Spring Constant of 0.1. 
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Figure 120. Plot of Ramp Exciting Force on all Three 
Pinned Connections at a Wave Angle of 180 Degrees 
with a Spring Constant of 0.1. 
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Figure 121. Plot of Ramp Exciting Force on all Three 
Pinned Connections at a Wave Angle of 225 Degrees 
with a Spring Constant of 0.1. 
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Figure 122. Plot of Ramp Exciting Force on all Three 
Pinned Connections at a Wave Angle of 270 Degrees 
with a Spring Constant of 0.1. 
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Figure 123. Plot of Ramp Exciting Force on the First 
Pinned Connection and Ship Heave Motion at a Wave 
Angle of 90 Degrees with a Spring Constant of 0.1. 
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Figure 124. Plot of Ramp Exciting Force on the First 
Pinned Connection and Ship Heave Motion Phase 
Angles at a Wave Angle of 90 Degrees with a Spring 
Constant of 0.1. 
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Plot of Barge Ramp Force Magnitude and Barge Heave Motion at 90 Degrees 
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Figure 125. Plot of Ramp Exciting Force on the 


Second Pinned Connection and Ship Heave Motion at 
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Figure 126. Plot of Ramp Exciting Force on the 
Second Pinned Connection and Ship Heave Motion 
Phase Angles at a Wave Angle of 90 Degrees with a 

Spring Constant of 0.1. 
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Plot of Barge Ramp Force Magnitude and Barge Heave Motion at 90 Degrees 
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Figure 127. Plot of Ramp Exciting Force on the Third 
Pinned Connection and Ship Heave Motion at a Wave 
Angle of 90 Degrees with a Spring Constant of 0.1. 
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Figure 128. Plot of Ramp Exciting Force on the Third 
Pinned Connection and Ship Heave Motion Phase 
Angles at a Wave Angle of 90 Degrees with a Spring 
Constant of 0.1. 
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Plot of Force on the ramp at 90 degrees, Vanous K Values 
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Figure 129. Plot of Ramp Exciting Force on the First 
Pinned Connection at a Wave Angle of 90 Degrees 
with Various Spring Constants. 
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Figure 130. Plot of Ramp Exciting Force on the 
Second Pinned Connection at a Wave Angle of 90 
Degrees with Various Spring Constants. 
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Figure 130a. Plot of Ramp Exciting Force on the 


Second Pinned Connection at a Wave Angle of 90 
Degrees with Various Spring Constants. 
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Figure 131. Plot of Ramp Exciting Force on the 


Third Pinned Connections at a Wave Angle of 90 
Degrees with Various Spring Constants. 
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Figure 131la. Plot of Ramp Exciting Force on the 


Third Pinned Connections at a Wave Angle of 90 
Degrees with Various Spring Constants. 
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V. CONCLUSIONS 


The main conclusions from this work are summarized below: 

1. |The comparison between a panel based three-dimensional hydrodynamic 
solver and a two-dimensional strip theory solver showed excellent agreement when 
modeling the LMSR. When the LMSR was modeled with the barges attached, although 
the strip theory solver was unable to account for the wave interaction affects, the added 
mass and damping coefficients as well as the wave exciting forces for the LMSR were 
still predicted fairly well. 

2. The use of the strip theory solver is more economical in that it requires 
significantly less computational time. Therefore, it may be effectively used to model the 
LMSR. A combination of two-dimensional and three-dimensional hydrodynamic solvers 
is the most efficient method to obtain the added mass and damping values for the 
LMSR/barge system. 

3. The strip theory solver was unable to accurately predict the hydrodynamic 
effects of the LMSR on the barge. This was attributed mainly to the irregular 
frequencies. 

4. Exciting forces that are induced due to proximity produce a pitching 
moment on the barges attached to the LMSR in beam seas. This moment causes the 
barges to pitch when they would naturally not exhibit any pitch motion. This results in 
larger barge motion than would be expected in beam seas and thus produces a large force 


on the pinned connections 1n beam seas. 
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5. Asheltering and magnification effect produces a dramatically lowered total 
force on the barges when seas are from 270 degrees and a significantly larger force when 
the wave direction is from 90 degrees. 

6. The king pin on the LMSR is the point where the most force is exhibited. 
This point is anywhere from 2 to 3 times larger than the connection on the barge. At 
higher frequencies, this force may be higher than the hydrodynamic forces which the ship 
experiences. 

7. With larger spring constants, the force acting on the pinned connections is 


significantly greater. 
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VI. RECOMMENDATIONS 


Based upon the data and conclusions, the following recommendations are made: 

1. Reformulate the problem as a four degree of freedom system. This would 
require modeling the ramp as non-nigid body. Preliminary studies indicate that the first 
torsional mode of the ramp is the most significant, so this would need to be incorporated 
first in the analysis. 

2. Utilize the transfer functions that were developed in this work to select the 
proper isolator properties. The minimization index would be the maximum stress in the 
ramp. 

3. Analyze all six motions of the LMSR and barges to verify that other motions 


are negligible as assumed in this problem. 
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APPENDIX A 
FIGURE OF LMSR/BARGE CONNECTION: 





Figure A. Basic Drawing of LMSR and Roll-On, Roll-Off Discharge Facility Connected 
together by Ramp. Black Dots Indicate the Locations and Number of the Pinned 
Connections. 
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